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We present a new way to compute thermodynamical observables on the lattice. We compute excited states 
and thermodynamical functions in the scalar model via the Monte Carlo Hamiltonian technique. We find agree- 
ment with standard Lagrangian lattice calculations, but observe lesser fluctuations in the results from the MC 
Hamiltonian. 



1. MONTE CARLO HAMILTONIAN 



The idea of the Monte Carlo Hamiltonian has 
been introduced in Refs.|Q-||]. Here we want to 
recall the basic idea, outline some of its char- 
acteristic features and present some selected re- 
sults. Comparing the Lagrangian approach with 
the Hamiltonian approach in lattice field theory, 
the Hamiltonian approach has not flourished in 
producing a wealth of numerical data like the La- 
grangian approach. On the other hand, there are 
a few topics, where the Lagrangian approach has 
only given slow progress, e.g. on excited states, 
wave functions, particle scattering. In our opin- 
ion, the lattice Hamiltonian approach might have 
done a lot better, if it were not for the lack of ef- 
fective many-body methods to solve such Hamil- 
tonian and compute physical observables. Recall 
that the Kogut-Susskind Hamiltonian of lattice 
QED3+1 and QC-D3+1 is a quantum mechanical 
operator in a high-dimensional Hilbert space and 
a numerical many-body solution is quite difficult. 

The Monte Carlo Hamiltonian using a stochas- 
tic basis is an attempt to address this problem. 
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Figure 1. Energy spectrum in a low energy win- 
dow. 



Let us outline the idea. First, the Monte Carlo 
Hamiltonian is not an operator of canonical form. 
It is rather defined in a subspace of Hilbert space 
and is given in terms of matrix elements. The 
physical idea is that this object plays the role 
of an effective Hamiltonian. It has a window of 
validity. This can be visualized as a window in 
the energy spectrum. Beyond this window, the 
effective Hamiltonian is unphysical. When com- 
puting thermodynamical functions, the window 
of the energy spectrum translates into a tern- 
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perature window. Such a window can be ob- 
served in the numerical data (see below). The 
notion of a window is familar in physics from 
the concept a scaling window and also reminds of 
Wilson-Kadanoff's renormalisation group, where 
a Hamiltonian is constructed which gives the cor- 
rect physics of critical phenomena (at phase tran- 
sitions), but is unphysical otherwise. 



trick is to write the path integral as a ratio). In 
order to get a reasonable representative subspace 
of Hilbert space (which also influences the size 
of the window), we took in quantum mechanics 
Xi and x / from a finite grid of position states and 
computed all possible combinations of transitions. 
Then an algebraic diagonalisation of such a ma- 
trix yields a finite set of eigenstates and cigcn en- 
ergies of the effective Hamiltonian, i.e. the Monte 
Carlo Hamiltonian. For test cases from quantum 
mechanics this has been shown to work |l| . 




Figure 2. Free energy F(f3). Monte Carlo Hamil- 
tonian (full dots) vs. standard Lagrangian lattice 

calculations (open dots). _ 

Figure 3. Same as Figj3, for average energy U((3). 



What is the guiding principle in constructing 
such effective Hamiltonian? In contrast to the 
transfer matrix formalism which defines a Hamil- 
ton operator in the limit at — > 0, we start by 
considering physical amplitudes corresponding to 
finite transition time T. The finiteness of T is 
in part responsable for the finite energy window. 
However, in contrast to Lagrangian lattice field 
theory, where one considers n-point functions cor- 
responding to vacuum-to-vacuum transitions, we 
go beyond the vacuum state and consider transi- 
tion amplitudes from other states which probe in 
part the wave functions of excited states. To give 
an example in quantum mechanics, we consider 
transition amplitudes G{xf,t = T;Xi,t = 0), 
where \xi > and \xf > denote eigen states of ini- 
tial and final position. We compute those quanti- 
ties via standard path integrals on the lattice us- 
ing Monte Carlo with importance sampling (the 



2. MANY-BODY SYSTEMS AND 
STOCHASTIC BASIS 

It is simple to convince oneself that a straight 
forward generalisation to many-body systems and 
field theory is not feasible. The solution is to 
select physically important degrees of freedom. 
Here we use as guiding principle the enormous 
success of Lagrangian lattice field theory, and in 
particular the successful working of Monte Carlo 
with importance sampling. In Lagrangian lattice 
field theory it is possible to estimate an observ- 
able by taking the mean of the observable over a 
relatively small number of "equilibrium' field con- 
figurations, which have been constructed using a 
weightfactor built from the action. 

We suggest that 'equilibrium' field configura- 
tions in Lagrangian lattice field theory have a 
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Figure 4. Same as Fig.||, for entropy S((3). Figure 5. Same as Fig.||, for specific heat C(/3). 



close analogy in Hamiltonian lattice field theory 
in terms of states of a stochastic basis. In other 
words, taking 'equilibrium' configurations at a 
fixed time slice yields Hilbert states (Bargman 
states) which represent physically important de- 
grees of freedom. This constitutes the so-called 
stochastic basis. Thus, via Monte Carlo with 
importance sampling, involving the weightfactor 
with the action, we generate a subspace of Hilbert 
space of a relatively small dimension (in the order 
of 100 to 1000), which represents the important 
degrees of freedom for the physics in a certain 
window. 

The demonstration that such repesentative ba- 
sis exists and allows to compute excited states 
and thermodynamical functions has been first 
made for the chain of coupled harmonic oscilla- 
tors (Klein-Gordon model). In Ref.|2||| the exci- 
tation spectrum has been compared with the an- 
alytic answer, establishing the existence of an en- 
ergy window and thermodynamic functions (spe- 
cific heat) have been computed. 

The method has been extended to the scalar 
model in Ref. |||. This model is not analytically 
soluble, hence we have chosen to compare our re- 
sults with standard Langian lattice calculations. 
Fig.|l] presents the the first 50 eigenvalues of the 
excitation spectrum of the Monte Carlo Hamilto- 
nian. A sudden increase of errors indicates the 
upper limit of the energy window. The thermo- 
dynamical functions free energy, average energy, 



entropy and specific heat are displayed in Figs.Q- 
^| The temperature window reaches from tem- 
perature zero to an upper temperature given by 
j3 ~ 1 — 1.5. The size of the window depends on 
transition time T, the size of the stochastic basis, 
statistical errors, etc. In Figs.|| - [s] one observes 
that the data from the Monte Carlo Hamiltonian 
show lesser fluctuations than those from the stan- 
dard Lagrangian approach. 
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